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Abstract 

Recent studies of the topological properties of a general class of lattice Dirac 
operators are reported. This is based on a specific algebraic realization of the 
Ginsparg- Wilson relation in the form 75(75!)) + (75-0)75 = 2a 2k+1 {^D) 2k+2 where 
k stands for a non-negative integer. The choice k = corresponds to the commonly 
discussed Ginsparg- Wilson relation and thus to the overlap operator. It is shown 
that local chiral anomaly and the instanton-related index of all these operators are 
identical. The locality of all these Dirac operators for vanishing gauge fields is 
proved on the basis of explicit construction, but the locality with dynamical gauge 
fields has not been established yet. We suggest that the Wilsonian effective action is 
essential to avoid infrared singularities encountered in general perturbative analyses. 



1 Introduction 

Recent developments in the treatment of fermions in lattice gauge theory are based on a 
hermitian lattice Dirac operator 75.D which satisfies the Ginsparg- Wilson relation[l]j 

l5 D + Dj 5 = 2aD l5 D (1.1) 

where the lattice spacing a is utilized to make a dimensional consideration transparent, 
and 75 is a hermitian chiral Dirac matrix. An explicit example of the operator satisfying 
(1.1) and free of species doubling has been given by Neuberger.0 The relation (1.1) led 
to an interesting analysis of the notion of index in lattice gauge theory. || This index 
theorem in turn led to a new form of chiral symmetry, and the chiral anomaly is obtained 
as a non-trivial Jacobian factor under this modified chiral transformation. Q This chiral 
Jacobian is regarded as a lattice generalization of the continuum path integral. The 
very detailed analyses of the lattice chiral Jacobian have been performed. || It is also 
possible to formulate the lattice index theorem in a manner analogous to the continuum 
index theorem. [0] ' [§ ' || An interesting chirality sum rule, which relates the number of zero 
modes to that of the heaviest states, has also been noticed. |[L0|| See Refs. [|TlH for reviews 
of these developments. 
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We have recently discussed the possible generalization of (1.1) and its implications. 
To be specific, we have discussed a generalization of the algebra (1.1) in the form 
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75(75^) + (75^)75 = 2a 2fc+i ( 75 ^) 2fe+i (1-2) 

where k stands for a non-negative integer and k = corresponds to the ordinary Ginsparg- 
Wilson relation. When one defines 



(1.2) is rewritten as 
or equivalent ly 
where we defined 



H = 7 5 aD (1.3) 

lb H + H lb = 2H 2k+2 (1.4) 
T 5 H + HT 5 = (1.5) 



T 5 = l5 -H 2k+1 . (1.6) 

Note that both of H and T 5 are hermitian operators. 

It has been shown that all the good topological properties of the overlap operator 



is retained in this generalization. |T2|, |13|| The practical applications of this generalization 
are not known at this moment. We however mention the charactristic properties of this 
generalization: The spectrum near the continuum configuration is closer to that of con- 
tinuum theory and the chiral symmetry breaking terms become more irrelevent in the 
continuum limit for k > 1. The operator however spreads over more lattice points for 
large k. 

2 Representation of the general algebra 

We first discuss a general representation of the algebraic relation (1.5). The relation (1.5) 
suggests that if 

H<f) n = a\ n <f) n , (<j) n ,<j) n ) = 1 (2.1) 
with a real eigenvalue a\ n for the hermitian operator H, then 

#(r 5 n ) = -aA n (r 5 n ). (2.2) 

Namely, the eigenvalues \ n and — A„ are always paired if A n 7^ and (r 5 n , r 5 0„) 7^ 0. 
We also note the relation, which is derived by sandwiching the relation (1.4) by <p n , 

(0«,750n) = («A n ) 2fc+1 for A n ^ 0. (2.3) 

Consequently 

|(aA n ) 2fc+1 | = |(0n,750n)| < 1 1 <\>n \ \ \ \^<\>n \ \ = 1- (2.4) 

Namely, all the possible eigenvalues are bounded by 

|A n | < -. (2.5) 
a 
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We thus evaluate the norm of r 5 0„ 

(r 5 n ,r 5 n ) = (0„,( 75 -^+ 1 )( 75 -^+ 1 )0 n ) 

= fa, (i - # 2fc+i 75 - 75 # 2fc+1 + H^y)<p n ) 

(2.6) 

where we used (2.3). By remembering that all the eigenvalues are real, we find that ra is 
a "highest" state 

T 5 0„ = (2.7) 

only if 

[1 - (a\ n ) 2 ] = (1 - aAn)(l + a\ n ) = (2.8) 

for the Euclidean positive definite inner product (0 n ,0n) = Y^ x 0n( x )0n(^)- 

We thus conclude that the states 4> n with A n = ±- are not paired by the operation T 5 cj) n 

and 

75-D0 n = ±-0n, 750n = ±0n (2-9) 

a 

respectively. These eigenvalues are in fact the maximum or minimum of the possible 
eigenvalues of H/a due to (2.5). 

As for the vanishing eigenvalues H<fr n = 0, we find from (1.4) that if 7 50 ra = 0, namely, 
H[(l ± 75 )/2]0 n = 0. We can thus choose 

7 5 £>0n = 0, 75 n = (j) n or 75 n = -<f> n . (2.10) 



To summarize the analyses so far, all the normalizable eigenstates <p n of 75-D = H/a are 
categorized into the following 3 classes: 

(i) n± ("zero modes"), 

l5 D<f) n = 0, 75 0„ = ±0n, (2.11) 

(ii) N± ("highest states"), 

75 .D0 n = ±-4> n , 750n = ±0™, respectively, (2.12) 
a 

(iii) "paired states" with < |A n | < 1/a, 

lbD<p n = \ n <p n , J 5 D(F 5 (j) n ) = -\ n (T 5 <f) n ). (2.13) 
Note that T 5 {T 5 (j) n ) oc (fi n for < |A n | < 1/a. 

We thus obtain the index relation 0. Q 
TrT 5 = ^(0 n ,r 5 n ) 

n 

= (#»> r 5^n) + i&n, T 5 0„) + Y (0»> T 5 0n) 

An=0 0<|A n |<l/a |A„|=l/a 



= J2 (&»> r 50«) 
A„=0 

= E (75 - ^ +1 )0n) 

A„=0 

A„=0 

= n + — n_ = index (2-14) 

where n± stand for the number of normalizable zero modes with 7s0 n = ±0 n in the 
classification (i) above. We here used the fact that T 5 n = for the "highest states" 
and that <p n and T 5 n are orthogonal to each other for < |A n | < \ja since they have 
eigenvalues with opposite signatures. 

On the other hand, the relation Tr7 5 = 0, which is expected to be valid in (finite) 
lattice theory, leads to ( by using (2.3)) 

7>7s = E(^>750n) 

n 

= E (0">750n) + Yl (<t>n,lb<t>n) 
A n =0 X„^0 

= n + - n_ + E (aAn) 2fc+1 = 0. (2.15) 

A„^0 

In the last line of this relation, all the states except for the "highest states" with A n = ±l/a 
cancel pairwise for \ n ^ 0. We thus obtain a chirality sum rule[10] 

n+ + N + =n- + N_ (2.16) 

where N± stand for the number of "highest states" with 750„ = ±0 n in the classification 
(ii) above. These relations show that the chirality asymmetry at vanishing eigenvalues is 
balanced by the chirality asymmetry at the largest eigenvalues with \X n \ — 1/a. It was 
argued in Ref.[14] that A^ ± states are the topological (instanton-related) excitations of the 
would-be species doublers. 

We have thus established that the representation of all the algebraic relations (1.2) 
has a similar structure. In the next section, we show that the index n + — n_ is identical 
to all these algebraic relations if the operator 75.D satisfies suitable conditions. 

3 Chiral Jacobian and the index relation 

The Euclidean path integral for a fermion is defined by 

J V^Vip exp[J ipDip] (3.1) 



where 



f^ = J2^)D(x,y)^(y) (3.2) 

x,y 
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and the summation runs over all the points on the lattice. The relation (1.5) is re-written 

as 

75 r 575J D + DT 5 = (3.3) 
and thus the Euclidean action is invariant under the global "crural" transformation [[| 

$(x) -> $'(x) = ${x) + i J24>(z)tl5F 5 (z, x)7 5 

z 

- il>'(v) = Hv) + * E ^(y, W W W ) ( 3 - 4 ) 

with an infinitesimal constant parameter e. Under this transformation, one obtains a 
Jacobian factor 

V^'Vijj' = JVipVip (3.5) 

with 

J = exp[— 2iTreT 5 ] = exp[— 2ie(n + — n_)] (3.6) 

where we used the index relation (2.14). 

We now relate this index appearing in the Jacobian to the Pontryagin index of the 



gauge field in a smooth continuum limit by following the procedure in Ref.|12[. We start 
with 

Tr{r 5 /(^£)} = Tr{r 5 /(^^)} = n + - n_ (3.7) 

Namely, the index is not modified by any regulator f(x) with /(0) = 1 and f(x) rapidly 
going to zero for x — > oo, as can be confirmed by using (2.14). This means that you can 
use any suitable f(x) in the evaluation of the index by taking advantage of this property. 
We then consider a local version of the index 

tr{T 5 f(^^)}(x,x) = tr{( l5 - H* k+1 )f((^)}(x,x) (3.8) 

where trace stands for Dirac and Yang-Mills indices; Tr in (3.7) includes a sum over 
the lattice points x. A local version of the index is not sensitive to the precise boundary 
condition , and one may take an infinite volume limit of the lattice in the above expression. 

We now examine the continuum limit a —>■ of the above local expression (3.8)[J We 
first observe that the term 

tr{H^f{^^)} (3.9) 

goes to zero in this limit. The large eigenvalues of H = a^D are truncated at the value 
~ aM by the regulator f(x) which rapidly goes to zero for large x. In other words, the 
global index of the operator TrH 2k+1 f(^j^) ~ 0(aM) 2k+1 -> for a -> with fixed 
M. 

We thus examine the small a limit of 

^{75/(^2)}- (3.10) 

lr rhis continuum limit corresponds to the so-called "naive" continuum limit in the context of lattice 
gauge theory. 
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The operator appearing in this expression is well regularized by the function f(x), and we 
evaluate the above trace by using the plane wave basis to extract an explicit gauge field 
dependence. We consider a square lattice where the momentum is defined in the Brillouin 
zone 

7T , 371 , 

<k u < — ■ 3.11 

2a ~ M 2a K J 

We assume that the operator D is free of species doubling, which is proved for the explicit 

construction of D\ in other words, the operator D blows up rapidly (~ -) for small a in 

the momentum region corresponding to species doublers. The contributions of doublers 

are eliminated by the regulator f(x) in the above expression, since 

trM(^)} ~ - (3.12) 

for a — > if one chooses f(x) = e~ x , for example. 

We thus examine the above trace in the momentum range of the physical species 

(313) 

We obtain the limiting a — ► expression 

limtr{7 5 /( ( - 7 ^ )}(x,x) 
= l imt r ^* c -**« 7te /((^!) c *to 

a^O /_£ (27r) 4 15J V M 2 ' 

= lim limtr f ^ e ^/((^) e ** 

,• j. f L d4ft -ikx ,/(^75-^)\ ikx 

= f{76/(^)} (3.14) 
where we first take the limit a — > with fixed fc M in — L < < L, and then take the 



limit L — > oo. This procedure is justified if the integral is well convergent. [12 1 We also 
assumed that the operator D satisfies the following relation in the limit a — > 

De ikx h{x) -> e ifcx (- ? - g fi)h{x) 

= i{ft + ig fi){e lkx h{x)) = ifl(e ikx h(x)) (3.15) 

for any /fcred (— ^ < k^ < | L ), and a sufficiently smooth function The function 

corresponds to the gauge potential in our case, which in turn means that the gauge 
potential A^(x) is assumed to vary very little over the distances of the elementary lattice 
spacing. 
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Our final expression (3.14) in the limit M — > oo reproduces the Pontryagin number in 
the continuum formulation (with e 1234 = 1)0 

Kmjr l5 f(^/M*) = tr^Cjh^l^f J ^J\-k^) 

= ■^trer#F ia ,F a p. (3.16) 

When one combines (3.7) and (3.16), one reproduces the Atiyah-Singer index theorem 
(in continuum R 4 space) .|7], || We note that a local version of the index (anomaly) is 
valid for Abelian theory also. The global index (3.7) as well as a local version of the index 
(3.8) are both independent of the regulator f(x) provided^] 

/(0) = 1, /(oo) = 0, f'(x)x\ x=0 = f'(x)x\ x=oo = 0. (3.17) 

We have thus established that the lattice index in (3.7) for any algebraic relation in 
(1.2) is related to the Pontryagin index in a smooth continuum limit as 

2 

n + -n^= f d 4 xJ-^tre^ af3 F^F af3 . (3.18) 

This shows that the instanton-related topological property is identical for all the algebraic 
relations in (1.2), and the Jacobian factor (3.6) in fact contains the correct chiral anomaly. 
(We are implicitly assuming that the index (3.7) does not change in the process of taking 
a continuum limit.) 

A detailed perturbative analysis of chiral anomaly for the general operators with k > 
has been performed, and the above result has been confirmed. \I2\ Also a numerical study 
of the index relation has been performed: The numerical result indicates the consistency 
of our analyses. [O 



4 Explicit construction of the lattice Dirac operator 

for k > 1 

We now comment on an explicit construction of the lattice Dirac operator which satisfies 
the generalized algebraic relation (1.2) with k > 0. We start with the conventional Wilson 
fermion operator Dw defined by 

D w (x, y) = ij^C^x, y) + B(x, y) - -m 5 Xty , 

(X 

C^{x,y) = ^[5 x+(i a,yU^(y) - 8 X! y +[la Ul(x)}, 

U„{y) = exp[ia(/A M (y)], (4.1) 
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where we added a constant mass term to D\y for later convenience. The parameter r 
stands for the Wilson parameter. Our matrix convention is that 7 M are anti-hermitian, 
(7**)t = — 7 M , and thus Q = 7 M C (U (n, m) is hermitian 

(? f (4.2) 

The Dirac operator for a general value of is constructed by rewriting (1.2) 
of relations 

H 2k+1 l5 + y 5 H 2k+1 = 2H 2 ( 2k+1 \ 

H 2 lb - lb H 2 = 0, (4.3) 

with H = ay^D. The second relation in (4.3) is shown by using the defining relation (1.4), 
and the first of these relations (4.3) becomes identical to the ordinary Ginsparg- Wilson 
relation (1.1) if one defines Hpk+i) = H 2k+1 . One can thus construct a solution to (4.3) 
by following the prescription used by Neuberger|Q] 

H ^ = h 5[1 + D w +1) i mm mr ] (4 - 4) 

where 

D§ k+1) = z((?) 2fc+1 + fi 2fc+1 - (I^+i (4.5) 

The operator H itself is then finally defined by (in the representation where H^k+i) is 
diagonal) 

H={H {2k+1) fl 2k+l (4.6) 

in such a manner that the second relation of (4.3) is satisfied. This condition (4.3) is shown 
to be satisfied in the representation where H^k+i) is diagonal. fT2|| Also the conditions 
< mo < 2r = 2 and 

2m 2k+1 = 1 (4.7) 

ensure the absence of species doublers and a proper normalization of the Dirac operator 
H. 



5 Locality properties of general operators 

We have explained that the general operators for any finite k give rise to correct chiral 
anomaly and index relations in the (naive) continuum limit. This suggests that those 
operators are local for sufficiently smooth background gauge field configurations. The 
locality of the standard overlap operator with k = has been established by Hernandez, 
Jansen and Luscher,|15[ and by Neuberger.||16|| 

As for the direct proof of locality of the operator D for general k, one can show it for 
the vanishing gauge field by using the explicit solution for the operator H in momentum 
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representation [O, ITS 



j fc + 1 ] fc + 1 / fc + 1 I fc S> 

H{ap,) = l5 (-)^(—=)^{(^Hl + M k )^-(^H^-M k )^TT^} 
1 a 



1 s fc+1 , 1 



fc+1 / ~ fc+1 



' 2 V^) 



75(tt) 2fc+1 (~^=) 2,1+1 {(V F (fc) + M k) 2k+1 - (\/F(k) - M fe )2*+i $ 

(5.1) 



where 



^ = (* 2 r +i +M fc 2 , 

M fc = E(l - c,)] 2k+l - m 2k+l (5.2) 



and 

s M = sin ap^ 
c M = cos ap M 

s/= 7 M sin ap M . (5.3) 

For = 0, this operator is reduced to Neuberger's overlap operator. [0 Here the inner 
product is defined to be s 2 > 0. This operator is shown to be free of species doublers for 
the parameter mo within the range < mo < 2 when we set r = 1, and 2m 2 l k+l = 1 gives 
a proper normalization of H, namely, for an infinitesimal p^, i.e., for \ap^\ <C 1, 

H ~ - l5 a p[l + 0{ap) 2 ) + 75 ( 75 a jf) 2k + 2 (5.4) 

to be consistent with H = ^^aD; the last term in the rigth-hand side is the leading term 
of chiral symmetry breaking terms. 

The locality of this explicit construction (5.1) is shown by studying the analytic prop- 
erties in the Brillouin zone.[12|] It is important to recognize that this operator is not 



ultra-local but exponentially local; fl7| the operator H(x, y) decays exponentially for large 



separation in coordinate representation 

H(x, y) ~ exp[— \x — y\/(2.5ka)). (5.5) 

An explicit analysis of the locality of the operator H^k+i) {not H itself )in the presence 
of gauge field, in particular, the locality domain for the gauge field strength | | has been 
performed. The locality domain for H-F^H becomes smaller for larger k, but a definite non- 
zero domain has been established. |12| The remaing task is to show the locality domain of 
H-F^H for the operator if = (H^ 2 k+i)Y^ 2k+1 ^ ■ Due to the operation of taking the (2/c+l)th 
root, an explicit analysis has not been performed yet, though a supporting argument has 
been given. jT|] . 
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6 Conclusion 



We have reported the recent investigation of topological properties of a general class of 
lattice Dirac operators defined by the algebraic relation (1.2). All these operators satisfy 
the index theorem and thus they are topologically proper. A precise proof of the locality of 
these general Dirac operators with fully dynamical gauge fields remains to be formulated. 
The operators with large k is expected to exhibit infrared singularities in perturbative 
analyses as is suggested by the construction of H^k+i) in (4.4), and thus the Wilsonian 
formulation of effective action, which is supposed to be free of infrared singularities, would 
be essential. 

Although we discussed only 4-dimensional theory, the recent developments in the 
treatment of lattice fermions[ll] may have some implications on 2-dimensional theory 
also, which is the main subject of this Symposium. In this respect, the fact that the 
lattice Dirac operators are not ultra- local but exponentially local fll5| may be of some in- 
terest. See Ref. [ffBJ for a Ginsparg- Wilson construction on a 2-dimensional fuzzy sphere. 
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